http://www.arocmag.com/article/02-2019-02-042.html 


高 速 铁路 移动 通信 下 基于 IMT-A 的 时 变 信 道 建 模 ” 
座 勇 "?， 攀 卓 宕 3， 沈 轩 帆 !， 胡 异 1， 姚 海 梅 !， 张 楠 1 


(1. 重庆 大 学 飞行 器 测控 与 通信 教育 部 重点 实验 室 , 重庆 400044; 2. 西安 电子 科技 大 学 综合 业务 网 理论 及 关键 技术 
国家 重点 实验 室 , 西安 710071; 3. 北京 大 学 信息 工程 学 院 , 广东 深圳 518055) 


摘 要 : 信道 建 模 及 其 仿真 方法 是 高 速 铁 路 移动 通信 系统 设计 的 难点 ， 同 时 已 有 的 面向 高 速 移动 下 IMT-A 时 变 信道 建 
模 的 研究 仍 存在 不 足 。 针 对 上 述 挑战 ， 提 出 了 一 种 高 铁 移动 下 具有 时 变 参 数 的 基于 IMT-A 信道 的 建 模 方法 。 首 先 ， 针 
对 IMTA Wo 区 场景 ， 设 计 了 一 个 马尔 可 夫 过 程 来 模拟 时 变 徐 的 数量 变化 ， 接 着 推导 出 随时 间 变 化 的 信道 参数 
的 表达 式 ， 其 中 包括 徐 的 延迟 、 功 率 、 离 开 角 和 到 达 角 ， 并 得 到 了 信道 冲 激 响应 ; 其 次 ， 分 析 了 信道 的 时 变 空 时 域 互 
相关 函数 、 时 变 自 相关 函数 以 及 平稳 间隔 ; 最 后 ， 对 所 提 信 道 模 型 的 统计 性 质 进行 了 仿真 ， 验 证 了 该 模型 具有 时 变性 
以 及 高 速 铁路 信道 具有 的 非 平稳 性 ， 并 且 证 明了 所 提 信 道 模型 可 用 来 模拟 高 速 铁路 信道 的 可 行 性 。 

关键 词 : 高 速 移动 ; 信道 建 模 ; 时 变 参 数 ; 马尔 可 夫 过 程 ; 非 平稳 性 

中 图 分 类 号 : TN918 doi: 10.3969/j.issn.1001-3695.2017.08.0887 


Time-varying IMT-A channel modeling for high-speed train mobile communication scenario 
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Abstract: Channel modeling and its simulation method is the difficult point of high-speed train mobile communication system 
design. At the same time, there are still some researches on IMT-A time-varying channel modeling for high-speed train mobile 
communication. Aiming at the above challenges, this paper presents a modeling method based on IMT-A channel with high- 
speed train mobile with time-varying parameters. First, for the urban macro-cell (UMA) scene of IMT-A, we design a Markov 
process to simulate the number of time-varying clusters, and then derive the expression of the channel parameters with time, 
including the delay, power, Angle and reach angle, and get the channel impulse response(CIR) . Secondly, the time-varying 
space-time cross-correlation function(CCF) , time-varying auto-correlation function(ACF) and stationary interval of channel are 
analyzed. Finally, the statistical properties of the proposed channel model are simulated, which proves the non-stationarity of 
the model and the non-stationarity of the high-speed train channel, and proves that the proposed channel model can be used to 
simulate the feasibility of high-speed train channel. 
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stationarity 
0 引言 媒体 通信 需求 。 众 所 周知 ， 信 道 是 通信 系统 的 基础 ， 它 的 模型 
一 及 仿真 方法 对 于 特定 场景 下 通信 系统 的 设计 具有 重要 作用 。 当 
近年 来 ， 随 着 高 速 铁路 的 发 展 ， 乘 客 对 高 质量 、 ee 前 面向 高 铁 场 景 的 标准 化 信道 模型 主要 有 IMT-advanced(IMT- 


数据 和 多 媒体 通信 提出 了 更 高 的 需求 。 然 ; A) 中 和 WINNER II 信道 模型 ， 它 们 均 假 设 信 道 满足 广义 稳定 
技术 还 无 法 满足 运行 速度 在 300 km/h 以 上 列车 的 高 质量 的 多 (wide-sense stationary, WSS) 并 且 忽略 了 其 非 平 稳 性 。 文 献 [2] 中 
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信道 测量 的 结果 表明 ， 标 准 信道 模型 明显 比 实际 测量 的 高 速 铁 
路 信道 具有 更 长 的 平稳 间隔 (定义 为 信道 满足 WSS 条 件 的 最 
大 持续 时 间 )， 而 WSS 假设 仅 对 实际 无 线 信 道中 的 短 时 间 间 隔 
有 效 。 综 上 所 述 ， 采 用 常规 平稳 信道 来 建 模 高 速 铁路 信道 的 方 
法 已 经 受到 日 益 严峻 的 挑战 。 
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的 已 条 路 径 中 的 一 条 和 0 条 子路 径 中 的 一 条 。 此 时 的 AoD 和 
第 一 单 射 / 禾 $ 有 关 ，AoA 来 自 最 后 的 单 射 / 禾 S$, 。 此外，a7 
和 or 分 别 是 基站 (Base Station, BS) 和 MS 的 视 距 (Line of Sight, 
LoS) 分 量 。 本 文 假设 BS 是 固定 的 ， 而 徐 5$ ， 以 w, 的 方向 和 速 
度 v 运动 。 当 本 文 考虑 对 于 每 个 路 径 发 生 多 次 散射 时 , 簇 9, 表 


前 ， 国 内 外 针对 高 铁 移动 通信 场景 下 的 信道 建 模 已 做 了 
量 研 究 ， 特 别 是 针对 非 几 何 随 机 模型 (non-geometrical 
stochastic model, NGSM)。 文 献 [3] 提 出 了 一 种 基于 高 铁 无 线 通 
信 的 有 限 状态 马尔 可 夫 链 的 NGSM， 使 用 马尔 可 夫 链 来 表征 移 
动 簇 (moving clusters, MC) 的 生 灭 性 质 。 所 提 模 型 能 够 利用 马尔 
可 夫 链 在 不 同 的 接收 信 噪 比 区 间 跟 踪 信 道 状 态 的 变化 ， 从 而 捕 
获 高 铁 无 线 信 道 随时 间 变 化 的 特性 。 但 是 ， 该 模型 的 非 平稳 的 
统计 性 质 很 难 通过 推导 得 出 。 此 外 , 文献 [外 提出 了 基于 NGSM 
的 非 平 稳 车 对 车 信道 。 基 于 在 高 架 桥 和 狭窄 通道 场景 下 的 高 铁 
信道 测量 数据 ， 文 献 [5] 提 出 了 一 种 有 限 状 态 的 马尔 可 夫 信 道 ， 
， 莱 斯 分 布 能 很 好 地 表征 小 尺度 衰落 在 两 种 高 铁 场 景 
下 的 测量 幅度 , 并 且 NGSM 可 有 效 地 捕捉 在 高 铁 信道 中 快 衰落 
的 动态 本 质 。 但 是 在 现 有 NGSM 的 研究 方法 中 ， 目 前 尚 没有 男 
向 IMT-A 时 变 信道 的 建 模 ,， 因 此 , 为 了 健全 和 完 
线 信道 建 模 理论 ， 研 究 这 类 符合 高 速 铁路 移动 通 
道 建 模 方法 具有 一 定 的 应 用 参考 价值 。 
为 此 , 本 文 提出 了 一 种 高 铁 移 动 通信 NGSM 场景 下 具有 时 
变 参数 的 基于 IMT-A 的 信道 建 模 方法 , 其 中 时 变 参数 包括 艇 的 
数量 、 功 率 和 延迟 ， 以 及 离开 角 (angle ofdeparture,AoD) 和 到 达 
角 (angles of arrival, AoA)。 首 先 本 文理 论 分 析 了 所 提 模 型 重要 
的 统计 特性 ， 包 括 局 部 空间 互相 关 函 数 (cross-correlation 
function, CCF)、 局 部 时 间 自 相关 函数 (auto-correlation function, 
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示 具 有 移动 方向 w。 和 速度 vy, 的 最 后 一 次 反射 。 对 于 MS， ays 
是 具有 速度 vy, 的 移动 方向 。 MC Ss ，MC 5s, 和 MS 的 移动 
分 别 使 用 向 量 功 ，V, 和 去 ,表示 在 图 中 。 此 外 , 假设 BS 和 MC 
5 之 间 的 初始 距离 D7, (4) 以 及 MS 和 MC 5s, 之 间 的 初始 距离 
DA () 是 已 知 的 。BS 和 MS 之 间 的 LoS 分 量 的 初始 距离 为 
) 。 信 道 参数 及 定义 通过 表 1 给 出 。 
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移动 端 ( 高 铁 列车 ) 


基于 IMT-A 信道 建 模 的 时 变 参 数 示 意图 


一 


于 MC 8 ，MC s, 和 MS We BS, MC 和 MS 
之 间 的 距离 将 随 着 时 间 上 变化 。 通 过 考虑 相对 速度 MC SS 和 
MC 5 的 概念 ， 本 文 可 以 假定 MS 只 相对 于 MC 5, 的 运动 。 


[ou 


ACF) 和 平稳 间隔 。 其 次 ， 通 过 MATLAB 仿真 验证 了 该 信道 模 
型 具有 非 平稳 性 ， 从 而 证 明了 理论 分 析 和 推导 的 正确 性 。 


1 ， 基于 IMT-A 的 信道 建 模 


在 高 速 移动 的 环境 下 ， 可 观察 到 的 MC 或 移动 端 (Mobile 
Station, MS) 的 快速 移动 ， 这 样 的 移动 变化 可 能 会 违背 前 面 所 提 
到 的 WSS 假设 。 在 这 种 情况 下 提出 具有 时 变 参 数 的 基于 IMT- 
A 信道 建 模 ， 几何 示意 图 如 图 1 所 示 ， 其 信道 冲 激 响应 (channel 
impulse respond, CIR) 可 以 表示 为 

h(t,7) = SOafe-z 人 (1) 


p=1 


et 


其 中 : 


h(D)= WD/Ox 


pn i sin(aoy (1)) +iOrk sin( on, (7) Ee 
g=1 


(2) 


这 jms|eos( py (D-Qus )ej9m 


于 MS 和 /或 MC 的 移动 ， 信 道 模 型 参数 P(1) ，z,(D)， 
酌 ,()，cz(D 和 cx (D 需要 使 用 适当 的 时 变 函 数 来 表示 。 其 中 


pa pa 


Qh) ， 02 (4)，AQ (1) 和 A (1) 分 别 表 示 初 始 AoD, 初始 


AoA，AoD 偏 移 和 AoA 偏 移 。 下 标 p，g 分 别 表示 宽带 信道 中 


从 信道 测量 获得 距离 D (1,)、Ds,(h) 和 Dios(4) 后 ， 则 时 间 
处 的 时 变 距 离 可 以 表示 为 式 (3)、(4) 和 (5)。 图 1 中 信道 参数 的 定 
义 如 表 1 所 示 。 
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表 1 图 1 中 信道 参数 的 定义 
Di,(?) BS 和 MC 5 之 间 的 距离 
Dw (人 MS 和 MC So 之 间 的 距离 
Dios (Dt) BS 和 MS 之 间 LoS 分 量 的 距离 
a 与 第 p 个 艇 /路 径 内 的 第 g 个 子路 径 相 关 的 


pa pa 


AoD 和 AoA 
LoS AoD 和 LoS AoA 


有 R 
Cros ， Cros 


Aew7 (1)， 


pq 


Aaw (1) 


AoD 偏 移 和 AoA 偏 移 


六 ， VD, 和 Ds MC S,，MC S, 和 MS 的 速度 矢量 


Vs A MC A 的 速度 和 MS 的 移动 方向 
Vp» QCp MC Z 的 速度 和 MS 的 移动 方向 
Vis’ Ons MS 的 速度 和 MS 的 移动 方向 
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Dss(?) = 
G3) 
JR ) + (VD —2. D3s(10) vt: cos(X + 0, (GO 一 cn) 


Dis (1D) 
， (4) 
(Das (10)) + —2: Dis (10) -vicos(os (10)— 0,) 
Dros(D = 
. ; a— (5) 
VD 和 (10) 十 (vist) -2 Dios (10) “Vst “COS(Qiss Qos) 
为 了 推导 出 所 需求 解 的 时 变 角 参 数 , 本 文 设置 四 个 辅助 角 ， 
即 图 中 的 pi (1) ，p ， 8 (1) 和 6 。 应 当 注 意 ， 在 图 中 ， 仅 展 


示 出 了 LoS 路 径 和 第 p 个 散射 传播 路 径 。 作 为 示例 , 艇 % 和 艇 
5; 分别 对 应 于 第 p 个 散射 路 径 的 第 一 和 最 后 的 跳动 艇 。 如 果 
MC s, 和 /或 MC 5s, 是 平稳 的 ， 则 MC 5s, 的 速度 yw 或 MC 5s， 
的 速度 可 以 被 设置 为 零 。 这 不 会 影响 该 基于 IMT-A 的 信道 
建 模 中 的 时 变 信 道 参数 的 生成 过 程 。 


其 于 上 述 描 述 , 可 以 通过 以 下 步骤 获得 基于 IMT-A 信道 建 
模 中 的 信道 参数 的 实现 ， 如 图 2 所 示 : 
生成 初始 信道 系数 
设置 汤 景 ”网络 布局 一。 分 本 传 所 条 件 。 | 一 让。 计算 路 径 吉 
生成 初始 乌 功 率 一 一 一 | 。 生成 初始 时 延 。 一 一 | 生成 相关 的 大 规 和 
生成 初始 信道 系数 一 一 | 绘制 随机 初始 相位 一 一 >| 进行 簇 的 随机 组 合 一 > 生成 初始 角 参 数 


生成 /更 新 信道 参数 


于 


生成 /更 新 簇 的 时 延 RE 一] 计算 新 生成 的 簇 ( 生 灭 过 程 》 FR 一 一 更 新 MS 和 MC 的 位 置 


不 


重复 每 个 村 


生成 /更 新 簇 的 功率 六 着 生成 /更 新 角 参 数 ”一 一 ”生成 信道 系数 ， 
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示 为 
ER=Ewc 二 EN (0) 
其 中 se = Roc (Flt he 0 
Bis = | (his Dar (8) 


其 中 : Rj 是 移动 艇 所 占 的 百分比 。 由 于 时 间 间 隔 较 短暂 ,每 个 
时 间 间 隔 w 一 t_ 可 以 被 假设 为 匀速 运动 人 1。 这 使 得 (7) 和 (8) 得 到 
进一步 简化 


Evyc =(t —ti) Ryec (vitv;) (9) 

Ens = (ti —ti 1 ) Vs (10) 

在 任何 4 时 刻 可 以 通过 在 1 时 刻 区 分 在 先前 CIR 中 已 经 

存在 的 新 生成 的 艇 和 存活 的 簇 。 通过 簇 的 生成 速率 p, 和 新 路 径 

的 复合 速率 p, 来 描述 这 个 马尔 可 夫 过 程 。 CIR 实现 中 的 簇 的 总 
数 的 期 望 也 被 定义 为 初始 数 ， 可 表示 为 


E{P()}=P(n)= 人 a (GD) 
ps 
观察 这 个 CIR 的 时 间 序列 ， 每 个 能 从 4 ,时 刻 的 一 个 CIR 
到 7 时 刻 的 随后 的 CIR 保持 概率 为 
R(e,) = er-Am (12) 


丸 此 ， 通 过 马尔 可 夫 过 程 产生 的 多 个 新 簇 的 期 望 为 


E{P 


new 


= 所 (1-epam 13 
| : ) 和 


经 过 演化 后 , 两 个 簇 之 间 的 相关 性 通过 se, 量化 。 若 =, 的 值 
较 高 ， 则 会 导致 5 时 刻 的 祖先 簇 的 属性 与 其 在 1 时 刻 的 继承 
簇 之 间 的 相关 性 降低 ,在 NLoS 的 城市 宏 小 区 (Urban Macro-cell， 
UMA) 场 景 下 ， 从 原始 IMTA 信道 模型 获取 的 初始 复数 
P(t,)=18 ,并 且 簇 的 每 次 改变 只 考虑 在 相 邻 状态 的 值 之 间 改 变 。 
基于 上 述 情况 ,本 文 进行 马尔 可 夫 生 灭 过 程 的 MATLAB 仿真 ， 
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图 2 基于 IMT-A 信道 建 模 的 信道 参数 实现 流程 


2 ”时 变 信道 参数 分 析 


2.1 时 变 簇 的 生成 


真 数 据 为 : p =0.6/m， p, =0.03/m ， vys =60m/s, v=15m/s， 


Vv =5m/s ，Ric =0.3。 仿 真得 到 的 时 变 艇 的 数量 P(1) 如 图 3 所 


20 T 

在 高 速 铁路 移动 通信 场 境 下 ， 随 着 MS 的 移动 ，MS 和 BS Re ] 
之 间 的 传播 环境 将 改变 。 此 外 ， 还 有 一 些 移动 散射 体 在 环境 内 
的 移动 也 会 对 传播 环境 造成 影响 。 因 此 ， 在 行驶 过 程 中 ， 一 些 
路 征 会 消失 ， 同 时 一 些 新 的 路 径 会 出 现 。 文 献 [6] 提 出 了 生 灭 过 。 。 鳞 | 
程 的 移动 久 的 建 模 方 法 。 在 时 变 情况 下 ， 铸 只 考虑 在 特定 时 间 三“ ] 
段 上 存在 。 随 着 时 间 的 推移 ， 新 的 徐 出 现 即 出 生 ， 并 且 保持 一 至 17s 
定 的 时 间 即 存活 ， 然 后 最 终 消 失 即 死亡 。 应 当 注意 的 是 ， 时 变 | L 
筷 是 指 本 文 假设 此 高 铁 信道 是 一 个 多 径 信道 ， 但 与 传统 信道 不 | 
同 的 是 ， 这 个 多 径 信道 的 径 数 P 会 随 着 时 间 变 化 ， 并 且 其 状态 a 

0 2 4 8 10 12 


转移 连续 ， 这 就 符合 马尔 可 夫 过 程 的 连续 状态 转移 的 行为 。 所 
以 ， 基 于 簇生 成 和 重组 的 过 程 ，MC(sxc) 的 移动 以 及 在 刀 -+ 和 
友 时 刻 的 时 间 间 隔 中 WMC(ewc) 的 移动 引起 的 信道 波动 可 以 表 


6 
时 间 t(s) 


图 3 时 变 艇 的 数量 分 布 


录用 稿 


从 图 3 中 可 以 看 出 , 艇 的 初始 值 为 18， 锡 
间 而 发 生 改变 ， 并 且 一 共有 16、17、18、 
数值 。 
2.2 时 变 延 迟 

随 着 MC (包括 MC s 和 MC s,) 和 MS 的 移动 , 第 p 
个 艇 的 传播 距离 将 是 时 变 的 。 时 变 延 迟 r, (1) 的 计算 如 下 : 


浇 的 数量 会 随 着 时 
19、20 这 5 种 情况 的 


7,(t)= (ADs, +AD )/ c+r, (1)= 


(14) 
(D5; (1) -Ds (1)) /e+ (Ds (7)— D& (1)) /ctr, (7) 
其 中 ，c 是 空气 中 的 光速 ，(ADI, +ADS,) 是 由 MS 和 MC 的 移 


动 引 起 的 传播 距离 的 改变 。 ADx = Di, (7)-Dss() 是 BS 和 MC 
3 之 间 的 距离 差 ，ADA = D& (ID)-D&(O) 是 MC S 和 MS 之 
间 的 距离 差 ，z(D 是 MC 8S 和 MC 9 之 间 的 虚拟 链 路 四 的 延 
迟 ， 可 以 通过 在 上 = 天 时 刻 使 用 一 阶 滤波 算法 计算 得 到 中 


人 ) 至 el ae, (a ) 于 (1 加 el A (15) 


其 中 工 服从 (Du /ct ) 分 布 ，x 是 取决 于 虚拟 链 路 和 场景 
的 相干 性 的 参数 ， ru 是 最 大 延迟 并 且 可 以 从 文献 [1] 获 得 ， 例 
如 ，r，=1885ns。 
值得 注意 的 是 ，z (4,) 可 以 计算 为 
t, (10)=7,(1)- (Dss (1)+ADS, (1))/e (16) 


其 中 初始 延迟 z (4,) 从 指数 延迟 分 布 中 随机 抽取 ， 如 在 原始 
IMTA 信道 模型 中 解释 的 给 定 情形 外。 将 式 (3)~(5) 代 入 式 (14)， 
获得 在 不 同 瞬时 时 间 + 的 第 p 个 簇 的 时 变 延 迟 即 r, (7) 。 然 后 通 
过 减 去 最 小 延迟 来 对 延迟 进行 归 一 化 并 按照 降序 排序 。 其 中 ， 
一 化 延迟 将 用 于 时 变 信道 系数 生成 ( 见 式 (1))， 但 不 用 于 簇 的 
生成 。 

2.3” 簇 的 时 变 功 率 

在 时 间 的 第 n 个 簇 的 随机 平均 功率 表示 为 W,()， 它 可 
以 从 原始 IMT-A 信道 模型 给 出 的 参数 产生 帆 。 因 此 ， 基 于 上 再 
对 时 变 延 迟 的 计算 ， 那么 在 时 间 t 的 第 p 个 簇 的 随机 平均 功率 
可 以 表示 为 


二 


(1 一 mps )(z (0)) 


ps "cps 


Wi() =e “om xl10 1 2. 
其 中 :js 是 延迟 分 布 比例 因子 ，os 是 延迟 扩展 ，& 是 每 簇 阴 
影 项 (以 dB 为 单位 )， 它 们 都 可 以 从 文献 [中 的 表 A1-7 中 获 
得 。 令 所 有 艇 功率 的 总 功率 等 于 1, 簇 的 时 变 功率 W,(i) 可 以 表 


示 如 下 : 


()/D "ws (18) 
2.4 时 变 AoD 
时 变 AoD 的 时 变 函 数 wz CD) 可 表示 为 
a) 如 果 


-AX < A (ho )<0,—x<o, < ,lh ) Uo (0) + 和 CQ < 元 
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那么 
OND = Q(t0) -Po lt) (19) 
如 采 -x<ar (4)<0， Qi(w)<o4<oer(4)+x， 那 么 
Q(t) = (to) + Bos (i) (20) 
b) 如果 
0<on lt) < 元 -TS Sh) IA) tA < 元 ， 
那么 
0 (D= 0 (10) + pa, 人 (21) 
如 采 0<ar (0)<x， oh(W)<o4<or()+x， 那 么 
QW,(D) = 0 (0)— Py, (1) (22) 
其 中 : 
pn (= 


Dss (tu)+ wtcos(o do) 一 os) 


(Dss (1)) 刘 (v7) +2D’, (to) vtcos (Qs, (0) 三 on 


(23) 
arccos 


2.5 时 变 AoA 
a) 如 果 


-A <a%(h)<0, -za, Sar (t) JU Atas (和 ) 和 ,< 元 ， 


那么 
oN (D) =0, -Bar(t) tx C4) 
如 果 -z<sas(o<0，oo(oO<sw< 和 区 +oo() ， 那 么 
or (DD)=0,+PB' (+x (25) 
b) 如 果 
0<on (hh) <x, As 0, so (0) I A+ oN (th) < a, <A, 
那么 
as,(D)=0,+PBs(t)+x (26) 
如 采 0<as (4)<x， ol(t)-x<o,<or(h)， 那 么 
oar (DD=0%, -P(t)+x (27) 
其 中 
pu 人 = 
vi— Dys (1 )cos(&, = on (10)) (28) 


arccos 


(DE (1)) + (vw) 2D (ojweos(w — os (1)) 


应 当 注 意 ， 式 (24)~(27) 中 的 “+r ”是 校正 器 ， 以 确保 结果 


在 范围 [x,x)。 


3 ”统计 性 质 


3.1 ”时 变 空 时 域 互 相关 函数 

对 于 所 提出 的 基于 IMT-A 信道 建 模 系统 ， 空 间 时 变 CCF 
不 仅 取 决 于 BS 处 的 天 线 间隔 Ac, 和 MS 处 的 天 线 间隔 A6, ， 
还 取决 于 时 间 t+。 因此， 该 信道 建 模 的 空间 CCF 被 称 为 局 部 空 
间 CCF ， 表 示 为 中 


Ri, (ACT ,AD ) = 
(29) 


(1/0 )>elw,() i A sin(a A aaar sin (a CD) 
gqg=1 
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3.2 ”时 变 自 相关 函数 
考虑 到 马尔 可 夫 生 灭 过 程 , 艇 从 yAr 到 上 存活 的 概率 表示 


为 尺 二 erPesarmcotosj 。 因 此 ， 第 p 个 艇 的 非 平稳 信道 模型 
策 局 训 而 ACF 可 以 表示 为 
六 (2 AD) = Ef{n, (nn, (中 全 
2 G0) 
Rax(1C)> BE{W (DJX(PADJ7(ADZ(AD 
gq=1 
其 中 : 
区 (4A1) ei sin (of, (1))-sin (as, (1+A1)) (3 1a) 
7 人 Ah) 本 ei sin(ah, (D)-sin(al, (+4n)) (3 1b) 
Z(t, Ah) = -ieos (eh (Wm OO) 有 eit eos (eh +A -0 (ra) G3 1c) 


3.3 平稳 间隔 


可 以 使 用 平均 功率 延迟 分 布 (APDP) 计 算 平稳 间隔 , 表示 为 
于 (LAD)= > G2) 
N 


其 中 : h(t,r)= Dh (4)5(r 
布 的 数量 ， 是 第 个 下 降 的 时 间 。 两 个 APDP 之 间 的 相关 系 


z)， 六 是 要 平均 的 功率 延迟 分 


数 可 以 计算 为 
0 
然后 将 平稳 间隔 计算 为 
T, (1)=max {Arlo, wyr] G4) 
其 中 7' 是 相关 系数 的 给 定 阐 值 ,具体 值 将 在 后 面 的 仿真 中 给 出 。 
4 ”仿真 分 析 
接 下 来 对 所 提出 的 高 铁 移动 通信 场景 下 的 基于 IMT-A 信 
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局 部 时 间 CCF 的 绝对 值 
5 


图 4 


所 


时 间 t(s) 


提出 的 基于 IMT-A 信道 建 


LE 模 的 三 维 局 部 空间 CCF 绝对 值 


空间 CCF 的 绝对 值 


部 


从 图 中 


间距 的 增 大 而 下 降 , 但 是 在 不 同 的 时 刻 , 波动 的 幅 
总 的 来 看 ，CCF 峰值 整体 是 平缓 略 有 增加 的 趋势 ， 并 且 


道 建 模 的 统计 性 质 进行 MATLAB 仿真 分 析 。 

4.1 ”时 变 空 时 域 互 相关 函数 

在 这 里 ， 本 文 使 用 NLoS 的 UMA 场景 和 主要 仿真 参数 如 
下 : A (to) =random, Di (n) 100，w, =15°, vy =30m/s, as (to) 


长 幅度 越 小 ， 
起 的 。 


结合 


可 以 看 出 , 虽然 


图 5 不 同 瞬 时 的 局 部 空间 CCF 的 绝对 值 


局 部 空间 CCF 是 随 着 MS 处 的 天 线 


最 后 波动 趋向 了 


合 图 4 和 5 的 分 析 ， 


信道 模型 有 


平缓 。 


度 并 不 


以 上 正 是 由 于 非 平 


相同 。 
时 间 越 
稳 性 引 


表明 该 高 速 移动 下 基于 IMT-A 的 
很 明显 的 非 平稳 性 。 
4.2 ”时 变 自 相关 函数 

继续 使 用 NLoS 的 UMA 场景 进行 仿真 ， 以 下 参数 用 于 仿 


=random, DR,(t)=150m, v=20m/s ， w, =120"。 真 局 部 时 间 自 相关 函数 ACF: cz (4)=random, DY (t,)=100， 
图 4 展示 了 所 提出 的 基于 IMT-A 信道 建 模 的 三 维 (3D) 局 部 。 Qa=15%，v =30m/s， a (10)= Tandom, Dys(h)=70m, &, ~ 
空间 CCF 的 绝对 值 。 从 图 4 可 以 看 出 , 随 着 MS 处 的 天 线 间 距 140”，v, =5m/s， Qs =120%， vs =20m/s。 
A5; 的 逐渐 增 大 ， 局 部 空间 互相 关 函 数 CCF 呈现 下 降 趋势 ， 并 6 表示 NLoS UMA 场景 的 基于 IMTA 信道 建 模 的 3D 
且 是 伴随 着 有 一 定 周期 性 波动 的 下 降 形式 ， 波 动 的 幅度 也 随 着 ”局 部 时 间 自 相关 函数 的 绝对 值 。 类 似 于 局 部 空间 互相 关 函 数 ， 
MS 处 的 天 线 间隔 的 增 大 而 明显 减 小 ， 最 后 趋 近 于 平缓 。 这 是 于 时 变 AoD 和 AoA, 局 部 时 间 自 相关 函数 的 绝对 值 随时 间 
由 于 该 高 速 移动 下 的 信道 的 非 平稳 特性 引起 的 。 变化 。 随 着 时 间 差 的 逐渐 增 大 , 局 部 时 间 自 相关 函数 ACF 呈现 
图 5 表示 从 图 4 截取 的 5 个 不 同 瞬时 即 大 0s、 左 1s、 顾 2s、 ”下 降 趋 势 ， 并 且 有 一 定 周 期 性 波动 的 下 降 形式 ， 波 动 的 幅度 也 
三 3s、 人 三 4s 的 局 部 空间 互相 关 函 数 的 绝对 值 。 随 着 时 间 差 的 增 大 而 明显 减 小 ， 最 后 趋 近 于 平缓 。 


图 7 表示 从 图 6 截取 的 6 个 不 同 瞬时 


即 扩 0s、 帮 1S、 


{=2s、 
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大 3S、 夺 48 和 三 5s 的 局 部 时 间 自 相关 函数 ACF 的 绝对 值 。 本文 


可 以 看 出 ， 局 部 时 间 自 相关 函数 是 随 着 时 间 差 的 增 大 而 下 降 ， 


增 大 ， 


且 在 不 同 的 时 刻 ， 波 动 的 幅度 并 不 相同 。 并 且 有 随 着 时 间 的 
波动 有 变 大 的 趋势 。 在 时 间 差 为 0.25s 时 ,不同 的 时 刻 的 
ACEF 基本 都 急剧 下 降 了 80% 以 上 。 随 着 时 间 差 的 增 大 ， 


波动 趋 


向 于 平缓 。 从 图 6 和 7 同样 证 实 了 该 高 速 移动 下 基于 IMT-A 信 


道 模型 的 非 平稳 性 ， 仿 真 结果 再 次 匹配 了 分 析 结 果 。 


局 部 时 间 ACF 的 绝对 值 


时 间 t(s) 


图 6 基于 IMT-A 信道 建 模 的 三 维 局 部 时 间 ACF 绝对 值 


4.3 


平稳 间隔 


0 0.25 0.5 0.75 1.0 1.25 1.5 1.75 
时 间 差 ，At(s) 


不 同 瞬时 的 局 部 时 间 ACF 的 绝对 值 


图 8 


显示 了 该 基于 IMT-A 信道 建 模 中 速度 为 v=90m/s 和 速 


度 为 *=20m/s 时 的 平稳 间隔 的 经 验 补充 累积 分 布 函数 (CCDF)09 


的 结果 图 ， 其 他 的 仿真 参数 为 ，f =930 MHz，y' =0.8。 对 于 该 
信道 建 模 ，D?,()=100，v =0m/s， DA (1) =70m，w, =0"。 


中 
QQ —e -v=90m/s 
—8— v=20m/s| ] 


口 o 
~ wm 


间隔 > 横 轴 ) 


陷 


概率 ( 平 
a 


0.4r 


eee 

1 1 1 1 

0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 
时 间 间 隔 (s) 


1 


0.3 


图 8 速度 分 别 为 =90m/s 和 w=20m/s 的 CCDF 


从 图 


8 可 得 ， 速 度 为 v=90m/s 时 ， 本 文 所 提出 的 信道 建 模 


的 平稳 间隔 的 CCDF 概率 如 下 : 当 概 率 为 0.8 时 对 应 的 时 间 间 
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隔 为 0.009s ， 以 及 概率 为 0.6 时 对 应 的 时 间 间 隔 为 0.02s; 速度 
为 v=20m/s 时 ， 其 平稳 间隔 的 CCDF 的 情况 是 : 当 概 率 为 0.8 
时 对 应 的 时 间 间 隔 为 0.018s， 以 及 概率 为 0.6 时 对 应 的 时 间 间 
隔 为 0.06s。 通 过 以 上 两 者 数据 的 对 比 可 以 看 出 ， 速 度 越 大 时 ， 
相关 性 越 小 , 从 而 证 明了 所 提出 的 基于 IMT-A 的 信道 建 模 方法 
的 可 行 性 ， 符 合 客观 实际 情况 。 


5 


道 的 建 模 方法 ， 使 ) 
基于 此 信 


结束 语 


本 文 提出 了 一 种 高 速 铁路 移动 通信 场景 下 基于 IMT-A 信 
马尔 可 夫 过 程 建 模 时 变 艇 的 生 灭 过 程 ， 并 
模型 ， 推 导 和 分 析 了 其 重要 的 统计 特性 ， 即 局 部 空 


证 


tt 


间 CCF 和 
ACF 和 平稳 间隔 的 MATLAB 进行 了 仿 
特性 ， 验 证 了 本 文 所 提 信 道 建 模 方法 的 可 行 性 


部 时 间 ACF 及 平稳 间隔 。 对 该 信道 建 模 中 的 CCF、 
， 发 现 其 具有 非 平稳 


Il 


o 
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